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Abstract 



The n-tree approximation scheme, introduced in the context of random di- 
rected polymers, is here applied to the computation of the maximum Lya- 
punov exponent in a coupled map lattice. We discuss both an exact imple- 
mentation for small tree-depth n and a numerical implementation for larger 
ns. We find that the phase-transition predicted by the mean field approach 
shifts towards larger values of the coupling parameter when the depth n is 
increased. We conjecture that the transition eventually disappears. 
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I. INTRODUCTION 



The development of analytical techniques to determine Lyapunov exponents in extended 
systems is certainly an important issue in view of the relevant information provided by 
them. One cannot, in general, expect to find exact solutions, as the problem is already non- 
trivial in the case of low- dimensional systems. Since most of the results published in the 
literature about Lyapunov exponents follow from numerical simulations, the development 
of effective "perturbative" techniques is very welcome as they can also provide information 
about interactions and correlations that are otherwise undetectable. 

More numerical results are certainly available for coupled-map lattices (CML), since the 
discreteness of both space and time variables allow simpler and faster simulations. The 
commonest coupling scheme for a CML is 

< +i = m 

y\ = ex\_ 1 + {l-2e)x\ + ex\ Jrl (1) 
The corresponding evolution equation in the tangent space is 

£ +1 = fW (eg-* + (1 " 2e)tf + eCli) • (2) 

from which one can see that even the computation of the maximum Lyapunov exponent 
(MLE) in a CML requires the simultaneous consideration of several issues: (i) space-time 
correlations of the local multipliers m\ = f f (yf); (ii) sign fluctuations of the multipliers 
which induce partial cancellations in the dynamics of the perturbation £*; (iii) correlations 
in tangent space induced by the spatial coupling. 

The third of the above issues is definitely the first to be clarified as it arises already in the 
presence of positive (^-correlated multipliers. It is precisely this problem that we shall address 
in the present paper, trying to determine the MLE in the random matrix approximation, 
i.e. assuming that all multipliers are independent, identically distributed, random processes. 
This is the standard assumption made in the study of Anderson localization in disordered 
systems so that we can say that our investigation can be naturally extended to such a case. 

A first attempt to determine the effect of spatial coupling has been made in [|I|], where 
the authors performed a mean field anlaysis, exploiting the analogy with the free-energy 
computation in directed polymers growing in random media. There, it has been found that 
the spatial coupling induces a shift in the value of the MLE from the quenched average 
(A = (logmj)) in the absence of coupling, to the annealed average (A = log(mj)) above a 
critical coupling value. 

In this paper we apply the so-called "n-tree" approximation scheme fHH, to obtain 
more refined analytical estimates of the MLE and to test the convergence properties for 
increasing depth of the tree. The growing evidence that many features of CML dynamics 
are indeed present also in chains of oscillators and in partial differential equations suggests 
that techniques developed for CML can be extended to such systems. 

The paper is arranged as follows. In section II we recall the essential lines of the n- 
tree approximation in directed-polymer theory and reformulate the approach in the present 
context. Section III is devoted to the numerical implementation, while the small coupling 
limit is investigated in section IV. Finally, in section V we present some remarks about the 
problem of estimating the MLE in a coupled map lattice and recall the open problems. 
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II. THE METHOD 



In this section, we first recall the n-tree approximation in the context of directed poly- 
mers, with reference to a (l+l)d structure. The approach is then explicitely formulated for 
the determination of the Lyapunov exponent. 

Let us consider all directed walks in a square lattice composed of the following displace- 
ments, j — > j + s, where j = (i,n) represents a generic site, while s G {(—1, 1), (0, 1), (+1, 1)} 
(see, e.g., Fig. 1), and attribute a random energy e(j) to each site. The statistical problem 
amounts to computing the partition function Z L (0) of all L-step walks departing from the 
origin, 

z L (o) = J2 e M-PE w ) , (3) 

w 

where the sum runs over all 3 L paths, E w is the sum of the energies in the sites visited 
during the walk, and (3 is the inverse temperature. 
One can write a recursive relation for Z^, 

Z i+1 (0)=exp(-/3e(0))$> L (s) (4) 

s 

with the initial condition Z Q (0) = 1. The free energy per unit length is nothing but the 
exponential growth rate of Z with L, F = — limi_ +oc (log Zl)/(L{3), where (■) represents the 
average over independent disorder realizations. 

The main difficulty preventing an exact solution of the above problem is the correlation 
among the partition functions appearing in the r.h.s. of Eq. (f|). The n-tree approximation 
scheme consists in iterating n times such a recursive relation (which automatically accounts 
for all correlations up to n steps), 



Z L +n(0) = eX P 



Sl,S2,...,S„ 



Z L (S n ) , (5) 



where Sj = Em=o,«-i s m, and in neglecting the remaining correlations, i.e. among the terms 
in the r.h.s. of the above expression. 

The problem can be solved exactly by introducing the generating function 

H L (x) = (exp{-e-^Z L }} (6) 

which satisfies the recursive equation 

H L+n {x) = ^ ft H L (x-^\ogM^ , (7) 

where Mj is the contribution of all paths of length n arriving at site j (see Fig. 1). The initial 
condition H = (exp{— e~^ x }) is a front interpolating the two fixed points P (H = 0), and 
Pi (H = 1) of Eq. (|7|). Since Pi is unstable, while P is stable, the front moves to the right. 
It can be easily seen that the front velocity is nothing but the free energy of the polymer. 
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The velocity can be determined by approximating the forefront as Hl(x) ~ 1— exp(— jx). 
Substitution of this Ansatz into Eq. (|7|) reveals that the front moves with a velocity that 
depends on 7, 

I IrninJ I ^ Imin 



where 



and jmin is the 7- value where G n takes its minimum 

dGM 



c?7 



. (10) 



The value of 7 is implicitely determined by the initial condition: by expanding H (x) for large 
x, one realizes that 7 = /3. Therefore, if /3 < 7 m j„, the free energy coincides with the annealed 
average of the weights Mj, while different 7-averages are selected at lower temperatures. This 
implies the existence of a thermodynamic transition occurring at (3 C = jmin- 

The analogy between the evolution equation (|2|) in tangent space, and the recursive 
relation (f|) for the partition function suggests that the whole procedure can be extended to 
the computation of the MLE. In this new context, the time t plays the role of the polymer 
length L and m\ can be considered as a quenched noise in a 2d environment. The growth 
rate of Z (i.e. the free energy), finally becomes the Lyapunov exponent. 

The differences with respect to the previous cases are the presence of the anisotropy 
factors e and 1 — 2e and the absence of a temperature. The first one is, in principle, only 
a technical variation which leads, however, to a strong "degeneracy" as we will show in 
the following. The temperature, instead, can be removed by setting (3 = 1; the role of 
the relevant control parameter will be played by the coupling strength e parameter. With 
the above simple indications, we find that the n-tree approximation A n of the Lyapunov 
exponent is 

a J G ni^tmin) 'fmin ^ 1 (~\~\\ 

A "-\G n ( 7 = l) 7mm >l ' 1 j 

The function G n (7) has the form 

Gn(7) = ^logj£((£(l^ U • (12) 

As shown in Fig. 1, Wj's are the directed walks on the lattice arriving at site j after n steps 
(the depth of the tree), k(wj) is the number of steps not involving a change of position, and 
rrii(wj) is the multiplier in the ith time step of the path Wj. 

Eq. QTID suggests the possible existence a phase transition upon changing e: if 7 min > 1, 
the Lyapunov exponent is given by the annealed average of the multipliers. The main diffi- 
culty in the implementation of this approach is the computation of G n and its minimization. 



4 



In the limit n — > oo, the approximation becomes exact: an interesting question concerns 
the convergence to the asymptotic value. Unfortunately, as n increases, the expression of 
G n quickly becomes so complicated that it is practically impossible to handle the analyti- 
cal expression. For this reason, in the next section we shall address the question from the 
numerical point of view. 

III. NUMERICAL RESULTS 

It is not only true that the analytical expression of the MLE becomes very complicated 
as n increases, but also an "exact" numerical implementation is not an easy task. Indeed, 
the average of disorder implies the computation of several multiple integrals and, even in the 
simple case of a uniform distribution of multipliers, the presence of a power 7 in Eq. (|T2|)) 
makes the integrals immediately undoable. The only case we have found, where it is possible 
to combine an exact solution of the integrals with a powerful numerical analysis, is that of 
a dichotomic distribution. Indeed, a generic multiple integral over K variables becomes 
a sum over all 2 K combinations of the variables, each properly weighted according to the 
probability of the two possible values of the multiplier. From Fig. 1, one can see that the 
number K of integrals to be performed is already equal to 7 for for n = 3 and j = — 1, which 
in turn requires summing up 128 different terms. As a result, even in this simple case, it is 
not possible to go beyond n = 5. 

Since a global rescaling of the multipliers yields a trivial shift of the MLE, we can assume, 
without loss of generality, that mi(t) = {1,6}. Moreover, for the sake of simplicity, and in 
order to maximize the effect of the fluctuations (which are responsible for the deviation of 
the MLE from the single-map case) we have assumed that 1 and b have the same probability 
1/2. 

The results for n = 1, . . . , 5, b — 3 and small e values are reported in Fig. 2, where a slow 
convergence towards the asymptotic value (numerically determined by iterating a chain of 
1000 maps) can be observed. In the inset of the same figure, one can also notice that the 
critical e-value, above which the Lyapunov exponent corresponds to the annealed average, 
steadily increases, in agreement with the numerical results that do not give any evidence for 
the existence of the high-temperature phase. 

Although an exact implementation of the n-tree scheme is unfeasible already for n > 5, 
one can consider it as a numerical algorithm. Indeed, one can imagine to iterate n times 
a perturbation initially localized in the origin. The amplitude in the site j represents an 
instance of Mj and the average required in Eq. (|12[) can be computed by summing over 
independent realizations of the stochastic process. 

Moreover, one can notice that 7 plays a similar role to q in the standard multifractal 
analysis; the only difference is that here, besides a local 7-average, a linear average over 
different sites is also required. It is therefore important to understand how 7 min behaves for 
increasing n, i.e. to clarify whether the actual value of the MLE does arise from a specific 
7 value. 

The only drawback of the numerical implementation is the need of a sufficient statistics, 
a constraint that becomes increasingly important for larger n-values, since an accurate deter- 
mination of Mj strongly depends on unprobable large deviations as usual in a multifractal 
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analysis. Notwithstanding this limitation, it has been possible to arrive at n — 50, much 
beyond the limits for an exact implementation. 

The results for e = 0.01 (reported in Fig. 3) confirm that the approximate values of the 
MLE approach from above the asymptotic value Aoo (denoted by the horizontal line in the 
figure). A numerical investigation of the behaviour of (A n — A^) versus n suggests that 
the convergence is presumably slower than algebraic. The slow variation of A n with n is 
confirmed by the poor improvement obtained by introducing the refined estimates 

K = nAn - {n ~ 5)A - 5 , (13) 
5 

(see the squares in Fig. 3), where the choice of "5" is simply dictated by the spacing of 
the numerical results. Notice that this procedure is very effective when the main finite-size 
effect arises from some rapidly decaying initial deviation, as it is the case of the Lyapunov 
exponent computed with the standard orthonormalization procedure Q. 

Besides allowing to compute the MLE, the n-tree approximation yields an estimate of 
the optimal 7 mi „ value. By comparing the results for the various depths, one finds that 
'jmin slowly decreases. We conjecture that ■j m in eventually converges to 0. The conclusion 
is suggested by the analogy between Eq. (|T2|) and the behaviour of the maximal comoving 
Lyapunov exponent A c (i>) ||, which is defined as the growth rate in the site i — vt of a 
perturbation initially localized in the origin. In a system with left-right symmetry, X c (v) 
reaches its maximal value for v = 0, where it coincides with the MLE. Therefore, for n large 
enough, the dominant contribution to G n (7) is given by the growth rates around the origin, 
i.e. by their logarithmic average. Accordingly, we expect that 7 min will eventually approach 
0. 

Such a tendency is confirmed for all e-values that we have considered, even well inside 
the supposed high-temperature phase where A n = G n (l). One such example is illustrated 
in Fig. 4, where we report G n {^) for e = 0.05 and different depths n. Notice that all G n {^f)- 
curves take the same value for 7 = 1; this is a consequence of the very definition of G n (j): 
independently of n, G n (l) coincides with the annealed average. For relatively small values 
of n, the minimum of G n is attained for 7 m m > 1 and the correct estimate of the Lyapunov 
exponent is given by the annealed average. However, upon increasing n, ^min steadily 
decreases until 7 min < 1 (in the case illustrated in Fig. 4, this happens for 20 < n < 30). 
The conjectured convergence of 7 min to zero implies the eventual disappearance of the phase- 
transition. 

Further light onto the n-tree approximation scheme can be shed by comparing it with 
a similar, though entirely heuristic, approach. The structure of G n (j) requires iterating 
an initially localized perturbation: it is therefore quite natural to consider a different initial 
condition, uniformly spreaded over 2n sites, and to iterate it for n steps by assuming periodic 
boundary conditions. At variance with the scheme understood in Eq. (]T2]), all sites are now 
statistically equivalent, so that we can estimate the MLE directly from their growth rate (or, 
better, from the average growth rate, to reduce statistical fluctuations). The performance 
of this empirical approach can be judged from the results reported in Fig. 3, where one can 
see that the convergence is now from below. Moreover, while the direct estimates are worse 
than the corresponding values obtained from the n-tree scheme, the opposite is true for the 
improved estimates. We can interpret this results as an indication that the actual value of 
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the MLE in an extended system follows from a delicate balance of several processes. A more 
effective reduction of finite-size effects can be presumably accomplished by introducing a 
still different definition of finite-time finite-space Lyapunov exponent. Whether one such a 
definition exists which applies to generic models is not obvious at all. 



IV. COUPLING SENSITIVITY 

One case that is worth investigating with the aid of the n-tree approximation is the 
small-coupling limit. In particular, it is interesting to study the scaling behaviour of the 
MLE for decreasing e. This problem has been already considered in Ref. with the help of 
a mean field approach. Here, we discuss the improvements arising from the implementation 
of the n-tree scheme. 

In this section, e will be always so small that even in the lowest approximation (n = 1), 
7 m i„ < 1 and the MLE is given by the low-temperature expression. In the 1-tree approx- 
imation, the Lyapunov exponent is given by A(e) = Gi(j m i n (e)), where G\ follows from 
Eq. (|12p with n = 1, 

G l { 1 ) = -\og{mT) + - log Ul-2ey + 2e<\ . (14) 
7 7 I J 

For small e, the above expression simplifies to 

Gxin) = -log(m^) + -log(l + 2^) . (15) 

7 7 

Since, for e — > 0, one must recover the value of the uncoupled case, 7 m j n (e) must go to 
zero in such a way that also e lmin ^ — > 0. Accordingly, one can further expand Eq. (|i~4|), 
obtaining 

G 1 ( 7 ) = (logm) + ^7 + 2- (16) 
2 7 

where r 2 = ((logm) 2 ) — (logm) 2 is the variance of the local expansion rate. With some 
algebra, it is possible to show that the value of •y m im minimizing Gi is given by 

lmim (e) = 2^ , (17) 
a 



where a = | loge| diverges for e — » 0. On the one hand, formula (|PTD confirms the correctness 
of the Ansatz jmini^) and that e 1 = l/a can be really considered a small parameter. 
On the other hand, we notice that the convergence to zero is extremely slow. 

By substituting (|i~7|) into the expression for G% and retaining only the leading terms in 
e, we obtain 



A 1 (a) = A + r 2 ^ + -^ , (18) 
a a log a 

where A = (log m) is the Lyapunov exponent of the uncoupled problem, the second term is 
responsible for the leading correction and the third smaller contribution is reported for the 
sake of completeness. 
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The above result, already derived in |J, has been here recalled because it allows intro- 
ducing all the key steps that are necessary to deal with higher-order approximations. In the 
following, we illustrate the case n = 2 with some detail and mention the further adjustments 
expected for larger depths of the tree. The complete expression for G n ( 7 ) is already rather 
complicated for n = 2, 

G 2 ( 7 ) = ^ log {A( 7 )} + ^ log {2^/1 ( 7 ) + 2^(1 - 2e)Va(7) + /sfr)} , (19) 
where 



A (7) = (™ 7 > 

/ 2 ( 7 ) = ((m 1 + m 2 ) 7 ) 

/ 3 ( 7 ) = ((eV + (1 - 2e) 2 m 2 + e 2 m 3 y) 

while m, mi,m2,m3 represent the multipliers to be averaged. 

With the aid of the same approximations made for n — 1, the above expression simplifies 

to 

G 2 ( 7 ) = (logm) + ^ 7 + - (20) 
2 7 

which, in turn, gives the following expression for the Lyapunov exponent 

A 2 (a) = A + T 2 ^ + ^-^— . (21) 
a 2 a loe; a 



By comparing Eqs. fl2~I|) and (IB), one can see that they differ only in the last term of the 
r.h.s., which is smaller by a factor two in the 2-tree approximation. The same is true (apart 
from the coefficient of the third term) for larger values of n. Accordingly, this analysis seems 
to indicate that the MLE grows as log | loge|/ log£, independently of the depth n. These 
conclusions are indeed confirmed by the numerical implementation of the n-tree scheme for 
several depths. This can be noticed in Fig. 5, where it is reported the behaviour of 

S n = (A n -A )|log£| (22) 

versus logs. The slow growth of is consistent with the expected log | log £:| behaviour. 
However, direct numerical simulations (full dots in Fig. 5) suggest that the MLE grows as 
1/loge, with no doubly logarithmic correction. In our opinion, the apparent contradiction 
can be solved by noticing that the determination of the correct scaling behaviour requires, 
in principle, to take first the limit n — > oo (to estimate correctly the MLE) and then the 
limit e — > 0. In the above analysis, we have, instead, exchanged the two limits. Since the 
determination of G n ( 7 ) requires estimating an exponentially growing (with n) number of 
contributions, it is reasonable to conjecture that increasingly small e-values must be reached 
before the leading term really overtakes the others. It is therefore conceavable that, before 
this asymptotic regime sets in, a different scaling region appears which becomes wider and 
wider for increasing n. The several limits involved in this process (we must not forget the 
role of 7 ) make a rigorous confirmation of this conjecture a rather delicate matter. Here, we 
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limit ourselves to recall that in two coupled Id and 2d maps a purely 1/ logs behaviour has 
been proved to arise ||. However, we cannot exclude that in the present case the scaling 
e-region is so small that it is has not been reached by our numerical simulations; it will be 
possible to give a definite answer only by developing a more effective perturbative technique. 

Finally, let us notice that the naive approach outlined in the second part of the previous 
section falls short of identifying the leading l/|loge| dependence, predicting only linear 
corrections in e. We suspect that this failure is due to the lack of an infinite time limit in 
the corresponding definition of the Lyapunov exponent. 

V. CONCLUSIONS AND PERSPECTIVES 

In this paper we have implemented the n-tree approximation for the computation of 
the MLE. The results have revealed a slow convergence towards the asymptotic value. We 
believe that more than being an intrinsic limitation of the method, this is an indication 
of the hardness of the problem. This feeling is indeed confirmed by the analogy with the 
Kardar-Parisi-Zhang (KPZ) equation. As already shown in 0, the logarithm h = log£ of 
the perturbation approximately satisfies the KPZ equation 

ht = eh x x + e(h x ) 2 + f (t) , (23) 

where the subscripts denote derivatives with respect to either time (t) or space (x) variables 
(assumed now to be continuous) and £ (t) is a noise-like term corresponding to the logarithm 
of the local multiplier. Accordingly, the average value of £ is nothing but the single-map 
Lyapunov exponent, while the MLE of the lattice is obtained by adding the average of 
the nonlinear term. Such a correction can be easily estimated by recalling that in Id, the 
probability distribution of h is exactly the same as for the linear Edwards- Wilkinson model, 
obtained by neglecting the nonlinearity M. The latter is the distribution of the standard 
Brownian motion, i.e. the product of independent Gaussians for the spatial derivatives h x 
f|. Since the variance of each Gaussian is inversely proportional to the coefficient of the 
Laplacian, it is seen that s(h 2 x ) is independent of e. Therefore, the KPZ equation provides 
for any coupling strength the same annealed average value, i.e. it even fails to find the low- 
temperature phase. This result indicates that one must go beyond the KPZ approximation 
of the tangent dynamics (adding higher-order derivatives and further nonlinear terms) if the 
e-dependence of the MLE is to be recovered. 

The identification of effective schemes to fasten the convergence of finite-size estimates 
of the MLE remains an open problem. Possible routes to be explored in the future are 
represented by corrections to the KPZ equation, or by suitable modifications of the n-tree 
approximation. 

Acknowledgments We thank the Institute of Scientific Interchange (I.S.I.) in Torino 
Italy, where this work was started. 
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FIGURES 



FIG. 1. Schematic diagram of all paths starting from site and arriving at site j in 3 steps 

FIG. 2. n-tree approximations of the MLE for n = 1, . . . , 5 (from top to bottom) compared 
with the numerical results (dashed curve). These and all the other results in the paper have been 
obtained for a dichotomic distribution of the local multipliers (m = 1,3, with equal probabilities). 
An enlargement of the region around the supposed phase-transition is reported in the inset. 

FIG. 3. Results of the n-tree approximation for e = 0.01 (bullets) and refined estimates as from 
Eq. ( |I3| ) (squares). Triangles and diamonds refer to the results of the "naive" approximation de- 
scribed in the end of section III and to the corresponding refinements, respectively. The horizontal 
line represents the MLE as determined with the standard numerical procedure. 

FIG. 4. G n vs. 7 for n = 20 (dashed line), n = 30 (dot-dashed line) and n = 40 (full line) 
obtained from the numerical implementation of the n-tree approximation with e = 0.05. 

FIG. 5. S n (see Eq. (p2|)) vs. e for different values of n. The solid lines correspond to 
n = 1, . . . , 5 (from top to bottom). Squares denote the results of the n = 40 approximation, while 
bullets correspond to the outcome of direct numerical simulations. 
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